


#include "PoreGasSourceKernel.h"
registerMooseObject("TrilobitaApp", PoreGasSourceKernel);
template<>
InputParameters validParams<PoreGasSourceKernel>()
{
  InputParameters params = validParams<Kernel>();
  params.addClassDescription("Compute Pore Gas Source kernel ");
  return params;
}

PoreGasSourceKernel::PoreGasSourceKernel(const InputParameters & parameters) :
  Kernel(parameters),
  _drhos_dt(getMaterialProperty<Real>("drhos_dt")),
  _porosity(getMaterialProperty<Real>("porosity"))
{

}

Real PoreGasSourceKernel::computeQpResidual()
{
//	std::cout << _k[_qp] <<std::endl;
  Real r(0);
	  r = _drhos_dt[_qp] * _test[_i][_qp];
	  return r;
}

Real PoreGasSourceKernel::computeQpJacobian()
{
	Real jac(0);
	return jac;

}







